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I. SUMMARY OF RESEARCH PERFORMED IN THE THIRD QUARTER 


Research activities in the past quarter have been focused upon 

(1) further modeling studies of the lo sodium cloud for region A and Region C, 

(2) development and preliminary testing of our lo potassium cloud model, 

(3) continued evaluation of magnetospheric plasma data for Saturn to be used 
in our Titan hydrogen torus model, and (4) collection and initial evaluation 
of Lyman a profile data for use in our comet particle trajectory model. 


1. Spatial Morphology of the Sodium Cloud 


Region A 

The east-west intensity asymmetry of Region A, discovered and documented 
by Bergstralh et al. (1975 and 1977), was recently attributed by Smyth (1983) 
to the effects of both solar radiation pressure (for very low velocity escape 
of sodium from lo) and the plasma torus sink (for higher velocity escape of 
sodium from lo). The need, noted by Smyth (1983), to further clarify the role 
played by the lo plasma torus sink, has been initially addressed this past 
quarter. This has been accomplished by using our new lo sodium cloud model, 
which explicitly includes the oscillating plasma torus sink, to calculate the 
brightness of Region A as viewed through the rectangular viewing slit adopted 
by Bergstralh et al. (1975, 1977). The effects of the oscillating plasma 
torus sink, the presence of an east-west electric field in the magnetosphere 
(Ip and Goertz, 1983; Barbosa and Kivelson, 1983), and the impact of the 
offset dipole nature of the planetary magnetic field (Acuna, Behannon, and 
Connerney, 1983) have been preliminarily assessed and are briefly discussed 
here. More quantitative results will be presented in the final report. 

The oscillation of the plasma torus about the satellite plane introduces 
a modulation in the brightness of Region A because of the varying sodium 
lifetime at lo, but it does not introduce any east-west intensity asymmetry 
(if the plasma properties are independent of system III magnetic longitude). 

An east-west electric field in the magnetosphere has recently been proposed by 
Ip and Goertz (1983) and Barbosa and Kivelson (1983) to explain the dawn-dusk 
asymmetry in the torus ion emission brightness measured by the Voyager UVS 
instrument (Sandel and Broadf oot , 1982; Sheraansky and Sandel, 1982). The 
electric field essentially shifts the plasma torus to a new center that is 
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displaced toward the east by about 0.18 to 0.20 Jupiter radii. For a given 
radial displacement from Jupiter, this provides an overall east to west shift 
of the plasma torus properties by twice this amount and gives rise to a signi- 
ficant east-west asymmetry in the sodium sink because of the sharp temperature 
gradient inside of lo's orbit. This asymmetry in the sink provides an east- 
west intensity asymmetry in the brightness of Region A that is in agreement 
with the value observed by Bergstralh et al. (1975, 1977). Inclusion of the 
offset dipole magnetic field in the sodium model provides a modulation of the 
east-west intensity asymmetry about the mean value established by the east- 
west electric field. This modulation may account for some of the apparent 
scatter present in the intensity data of Bergstralh et al. (1975, 1977). 

Region C 

The directional features of the sodium cloud, discovered by Pilcher (see 
Hartline, 1980), extend from Region B into Region C. Within the past year, 
changes in the north-south directions of these features have been successfully 
understood through our model-data comparisons which have involved an active 
collaborative effort with Pilcher. The extremely interesting results of this 
collaborative effort are presented in a paper submitted in the past quarter to 
The Astrophysical Journal for publication. A copy of this paper is attached 
as an appendix to this progress report. Our paper is to be published jointly 
with a companion paper entitled "Ejection of Atoms and Molecules from lo by 
Plasma-ion Impact" submitted by E.M. Sieveka and R.E. Johnson of the Univer- 
sity of Virginia. 

In short, the directional features can result from a source of high velo- 
city (''ZO km sec~M sodium combined with the oscillating neutral sodium sink 
provided by the plasma torus. The phase relationship between the features' 
directions and lo's magnetic longitude can be understood if escaping sodium is 
initially directed at near right angles to lo's orbital motion. The initial 
directions and speeds of sodium atoms escaping lo to form the directional 
features can be understood in terms of a magnetospheric-wind-driven escape 
mechanism. The companion paper submitted by E.M. Sieveka and R.E. Johhnson 
addresses the subject of ejection of sodium atoms from lo and provides support 
for the magnetospheric-wind-driven escape mechanism. 
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2. lo Potassium Cloud Model 


The development of a model for the lo potassium cloud has been completed 
this past quarter. The model contains all the refinements present in our most 
recent lo sodium cloud model. The lifetime of potassium, determined for elec- 
tron impact ionization, may be specified for plasma conditions reflecting the 
encounter times of Voyager 1 and Voyager 2 with Jupiter. In calculating the 
brightness of the cloud as well as the acceleration effects of solar radiation 
pressure, the solar spectra for the 7665 A and 7699 A lines of potassium as 
determined for our model by Kurucz (1982) and reported in our last annual 
report (i.e., June 1, 1981 to May 31, 1982 period) were adopted in the model. 
The model has undergone initial testing and preliminary model runs are antici- 
pated in the next quarter. 

3. Titan Hydrogen Torus Modeling 

Effort this quarter has been limited to examining the two-dimensional 
(radial and vertical) density and temperature information for ions and elec- 
trons in Sturn's magnetosphere supplied recently by Sittler (1983). This 
information will be incorporated in the Titan torus model in the last quarter 
to calculate the lifetime of hydrogen in the circumplanetary environment. 

4. Modeling of Cometary Hydrogen 

The full solar disk Lyman-a profile as determined by Lemaire et al. 

(1978) has been adopted for use in our comet particle-trajectory model. This 
profile, shown in Figure 1, is the only available high-resolution full disk 
Lyman-a profile of the Sun (Skumanich, 1984). The profile has been put in 
digital form for use by the model. In our model calculations for the extended 
hydrogen cloud, the profile will be utilized to calculate the radiation pres- 
sure acceleration of H atoms and their brightness produced by solar resonance 
scattering. The process of incoporating the profile into the model, initiated 
in the past three months, will be completed during the next quarter. 
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II. Program for the Last Quarter 


Research activities in the last quarter will involve (1) implementing the 
Saturn magnetospheric plasma data in our Titan torus model and performing 
preliminary model calculations, (2) completing the incorporation of the Lyman- 
a profile in the comet particle trajectory model and performing initial model 
calculations for the extended H atmospheres of comets, and (3) performing 
preliminary model calculations for the lo potassium cloud. 
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ABSTRACT 


Elongated features in lo's sodium cloud, directed away from Jupiter and 
both to the north and south of the satellite's orbital plane, have been ob- 
served. The north/south directions of the features are correlated with lo's 
magnetic longitude, suggesting a formation mechamism involving the oscillating 
plasma torus. We show by means of a model analysis that the features can re- 
sult from a source of high velocity (~ 20 km sec”^) sodium combined with the 
oscillating neutral sodium sink provided by the plasma. The phase relationship 
between the features' directions and lo's magnetic longitude can be understood 
if escaping sodium is initially directed at near right angles to lo's orbital 
motion. The directionality of the features requires that the sodium flux from 
equatorial regions be higher than that from the poles. The initial directions 
and speeds of sodium atoms escaping lo to form the directional features can be 
understood in terms of a magnetospheric-wind-driven escape mechanism. The one 
sequence of directional feature observations that we have analyzed in detail 
implies a high speed sodium source rate of ~ 10 ° atoms sec . 
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I . INTRODUCTION 


Sodium D-line emission in the spectrum of lo, Jupiter's innermost 
Galilean satellite, was discovered by Brown (1974; Brown and Chaffee, 1974). 
Shortly thereafter, it was determined that this emission emanates from a sub- 
stantial volume of space surrounding the satellite (Trafton et al . 1974; 

Mekler and Eviatar, 1974; Wehinger and Wyckoff, 1974). The morphology of the 
sodium cloud, discussed in the review by Pilcher and Strobel (1982), leaves no 
doubt that the sodium originates on lo. The directions and velocity disper- 
sion of the atoms ejected from lo to populate the cloud have been a topic of 
considerable study, involving the analysis of both sodium D-line profile data 
and spatial intensity data. No escape mechanism has yet been identified, how- 
ever, that can properly explain both these spectral and spatial signatures. 

Trafton (1975) and Trafton and Macy (1977) have shown that the center of 
the sodium D-line, measured with lo in the field of view, exhibits a wave- 
length shift which depends upon the orbital phase of the satellite. In addi- 
tion, the line profile shows an asymmetric wing (or skirt) extending to dop- 
pler shifts relative to lo of up to 18 km sec~^. The skirt appears to be 
stronger when lo is near the magnetic equator. The asymmetric skirt is seen 
on the blue side of the line center when lo is near eastern elongation and on 
the red side when lo is near western elongation, so that the sodium atoms pro- 
ducing it always appear to be traveling more rapidly in the direction of lo's 
motion than the satellite. In analyzing these observations, Trafton and Macy 
(1978) concluded that the asymmetric skirt was produced by high-speed sodium 
( <20 km sec”^) streaming from lo, preferentially in the direction of the 
satellite's orbital motion, and that the line core could be produced either by 
an atmosphere stationary with respect to lo or by a cloud expanding primarily at 
low velocities (<3 km sec"^). They showed that a speed dispersion composed of a 
primary low-velocity peak for the line center and a secondary (but broader) 
high-velocity peak for the asymmetric skirt was sufficient to reproduce the 
characteristics of their data. 

Smyth and McElroy (1977) presented line profile calculations for lo at 
eastern elongation to explain the general characteristics observed in the data 
of Trafton (1975). They concluded that sodium ejected primarily from the lead- 
ing inner quadrant of the satellite at velocities near 3 km sec ^ would explain 
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the shift of the line center and that a speed dispersion ranging up to 15 km 
sec~^ would produce the asymmetric skirt. Carlson et al. (1978) analyzed three 
profiles: one line profile which they measured 8 arc seconds from lo with the 

satellite almost 45 degrees past western elongation, and two different line 
profiles measured by Trafton (1975), one of which was independently analyzed by 
Trafton and Macy (1978). Carlson et al. concluded that these three observa- 
tions were best matched by sodium atoms ejected from lo's leading hemisphere 
with a velocity dispersion that must be similar to the sputtering distribution 
that they adopted. The velocity distribution they assumed peaks at ~3 km sec ^ 
and extends to velocities as high as those inferred by Trafton and Macy (1977). 
Using the model of Carlson et al . (1978), Porco, Trauger and Carlson (Trauger, 
1980) also analyzed a line profile that they measured near the satellite when 
lo was at western elongation. They found that their observation was explained 
by radial sodium ejection from the satellite hemisphere centered on the sub- 
Jupiter point (i.e., the inner hemisphere). 

Although the exact angular and velocity dispersions of sodium are not very 
well defined by the rather broad hemispherical source regions adopted in the 
above analyses, certain primary conclusions emerge. To explain the line pro- 
file shapes, sodium ejected from lo must be described by a velocity distribu- 
tion that includes speeds up to 20 km sec~^; many of these high-speed atoms 
must have a significant velocity component in the direction of lo's motion. 
Since the preferential direction of these velocity components is in the same 
direction as that of the magnetospheric plasma moving past lo, and since the 
size of the asymmetric line profile wing appears to be correlated with lo's 
magnetic latitude, a dependence of the sodium source on the magnetospheric flow 
is suggested. Furthermore, from the direction and magnitude of the high-speed 
sodium, one would also expect to observe sodium at a considerable distance from 
lo. Sodium emission has in fact been reported in Jupiter's equatorial plane at 
distances as large as ~35 Rj (Pilcher and Schempp, 1979) and ~53 Rj (Trafton 
and Macy, 1978) from the planet. These distances should be compared to lo's 
orbital radius of ~6 Rj. 

More recent D-line profile data have been obtained by Brown and Schneider 
(1981) at distances up to 20 Rj from Jupiter. In most spectra they observed a 
narrow emission feature which they concluded was due to sodium atoms on bound 
orbits with apojoves in the field of view and perijoves near lo's orbit. Ejec- 
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tion from lo in the forward direction at velocities relative to the satellite 
of 2-4 km/sec would yield atoms on the inferred orbits. On two occasions, how- 
ever , Brown and Schneider observed a broad emission feature whose central wave- 
length corresponded to a velocity of ~30 kra/sec relative to Jupiter and which 
exhibited a tail extending to atomic velocities of <100 km/sec in Jupiter's 
rest frame. They attributed this feature to some combination of charge ex- 
change between neutral and corotating ionized sodium, and elastic collisions 
between sodium atoms and corotating heavy ions. They noted that the rates of 
both of these processes should be at a maximum when To is near the symmetry 
plane of the plasma torus, and that their two observations projected back in 
time to this geometry. Brown and Schneider noted that although the velocities 
they measured were considerably higher than those reported by Trafton and Macy, 
their observations and those of the latter authors were likely related. 

Spatial intensity measurements of the sodium region-B cloud in lo's imme- 
diate vicinity appear to provide information primarily on material that left lo 
at relatively low velocities. Matson et al. (1978) presented two images show- 
ing a banana-shaped cloud surrounding lo with more sodium atoms preceding the 
satellite in its orbit than following. They found they could account for this 
shape by assuming sodium atoms were ejected with a velocity distribution 
characteristic of ion sputtering from a hemisphere intermediate between lo's 
leading hemisphere and that facing Jupiter. Smyth and McElroy (1978), analyz- 
ing a set of 56 images of the region-B cloud obtained by Murcray and Goody 
(1978; also see Murcray, 1978), concluded that sodium appears to originate on 
the inner hemisphere or some combination of the inner and trailing hemispheres 
of lo. They found that the data could be matched by models with a relatively 
constant hemispheric source of 2 x 10 atoms sec , ejected at a characteris- 
tic velocity of 2.6 km sec”^ relative to lo. Their results indicate that the 
data of Matson et al. (1978) can probably be accounted for by the low velocity 
peak in the sputtering velocity distribution. An improved lo sodium cloud 
model including the effects of the oscillating plasma torus, used in this paper 
to analyze the newly discovered sodium directional features, has recently been 
applied to the region-B cloud (Smyth and Combi, 1983). Because of the inclu- 
sion of this severe, spatially non-uniform sink for sodium, it now appears that 
the observations of Murcray and Goody (1978) can be explained by low-velocity 

sodium ejected symmetrically from lo's inner and outer hemispheres with an in- 

—1 

creased source rate of about 2 x 10 ° sec 
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In this paper we report the results of an imaging study of the sodium 
D 2 ~line intensity in the vicinity of lo which shows the effects of high velo- 
city ejection from the satellite. The manifestation in our data is in the form 
of previously unrecognized directional features. We show that the features, 
which are directed away from Jupiter and appear to oscillate north and south of 
lo's orbital plane, may be explained in terms of a high velocity sodium source 
and an oscillating sink provided by the lo plasma torus. The directions and 
velocities of sodium atoms required to account for these features provide 
direct evidence for a magnetospheric-wind-driven gas escape mechanism from 
lo. We also show that during a portion of February 1980, the sodium emission 
intensity surrounding lo, including that of the directional features was 
anomalously high. The difference with prior and subsequent observations 
indicates that this was a period of unusual activity on or near the satellite. 
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II. OBSERVATIONAL AND DATA REDUCTION TECHNIQUES 


We obtained sodium D 2 ~line images using a 2 .5- A bandwidth (FWHM) inter- 
ference filter and an image-tube/photographic-plate detector on the 2.2-m tele- 
scope at Mauna Kea Observatory in the manner described by Pilcher (1980). A 
14-arcsecond diameter occulting disk in the form of 1.5-mm aluminum dot de- 
posited on a window in the telescope focal plane was used to prevent light from 
lo's disk from entering the optical system. The field of view was 110 arc- 
seconds in diameter, which corresponded to 5.0 to 5.5 Rj depending upon the 
geocentric distance of Jupiter at the time of the observation. Exposure times 
were generally three minutes. The interference filter was designed so that at 
0°C and normal incidence its passband was centered slightly (~4 A) to the red 
of the rest wavelength of the sodium D 2 line (5890.0 A). We tilted the filter 
to tune it toward the blue for the observations. 

It is possible to do these observations with or without a collimator pre- 
ceding the filter in the optical train. If a collimator is not used, then the 
divergence angle of the beam emanating from any image point in the telescope 
focal plane will cause different portions of the beam to be transmitted by the 
filter at somewhat different wavelengths. In this application the result would 
be a loss in contrast between the monochromatic image and the background con- 
tinuum. Since we wished to detect the faintest radiation possible, we chose to 
collimate the beam from the telescope. This, however, introduced another vari- 
ation. Collimated bundles from different points in the telescope focal plane 
do not strike the filter at the same angle. For a monchromatic source, the 
effect is a variation in the transmission across the field of view. This is 
illustrated in Figure 1, in which we show a sequence of three-minute exposures 
taken at filter tilts separated by 1°. These images were taken when lo was 
within a few degrees of western elongation. In Figure 2 we show contour plots 
of the system spatial transmission corresponding to three of these images. The 
spatial transmission functions were derived from room temperature, laboratory 
tracings of the filter transmission and an empirical correction for the change 
in central wavelength with temperature determined from a comparison of the lab- 
oratory and Jovian data. It is clear from Figures 1 and 2 that although the 
transmission in this case at filter tilt angles of 1° or 2° is sufficiently 
uniform to show sodium emission over the entire field, this is not so at higher 
tilt angles. Since approximate spatial transmission functions are known, it is 
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possible to correct images for this effect in all areas where sodium emission 
has been detected and the system transmission function is significantly above 
zero. All of the contour plots presented below have been corrected in this 
manner . 

The data were traced on a Perkin-Elmer PDS lOlOA microdensitometer con- 
trolled by means of a DEC PDP 11/40 minicomputer. The digitization interval 
was ~0.5 arcseconds and the scanning aperture width was slightly less than the 
spacing between samples. For later processing and display in the form of con- 
tour maps, 4x4 pixels of each image were averaged. This filtering yielded 
higher signal-to-noise ratio, but a resolution of only ~2 arcseconds or 
~0.1 Rj. Since the seeing and telescope tracking were generally better than 
this value, this is the approximate spatial resolution of the digital data in 
final form. 

Before filtering or background subtraction (see below) we converted ap- 
proximately from photographic density to relative intensity using calibration 
curves obtained during spectroscopic observations on other nights. Calibration 
curves obtained on several nights did not vary significantly from one another 
at relative intensities below 400 (arbitrary units). At higher relative inten- 
sities, however, the curves sometimes diverged, and in no cases did the curves 
extend beyond 600. Our relative intensity scale (uncorrected for system spa- 
tial transmission) should therefore be regarded as only an estimate at values 
much above 400. This statement applies particularly to data obtained during 
February 1980 in which values exceeding 1000-2000 are common. For images such 
as these, in which the maximum photographic density exceeds the maximum density 
on the calibration plates, we extrapolated the calibration curves. The inac- 
curacies inherent in this procedure render these high relative intensity values 
useful only as indicators of extremely bright emission. 

In general, we obtained images in sequences of increasing filter tilt 
angle. The angular increment was 1° and the sequences usually extended to 
angles (~6°) at which essentially no D-line radiation was transmitted. The 
best such image for each night was used as a measure of the background and sub- 
tracted from each D-llne image obtained on that night. All of the contour 
plots presented here are of background-subtracted data. 

The background images obtained to the east of Jupiter exhibit an essen- 
tially monotonic decrease in intensity with increasing distance from the 
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planet. Those obtained to the west, however, exhibit a sudden increase in the 
background intensity across a line on the Jupiter side of the field roughly 
parallel to the projected Jovian rotation axis (cf. Fig. If). This discontinu- 
ity appears to be the result of an internal reflection from the telescope ter- 
tiary mirror which in its stowed position introduces an east-west asymmetry 
into the scattering characteristics of the telescope. 

Another characteristic of background images obtained on both sides of the 
planet is lower intensity at the southern edge of the field than at the north. 
This does not appear to be due to the flat field response of the image tube, 
which is known not to have variations in the red exceeding 10% (Stockton, pri- 
vate communication). There are two possibilities. One is that the north-south 
variation results from non-uniform scattering of the Jovian background radia- 
tion. In this case the effect on the sodium images will be completely removed 
in the process of background subtraction. The second possibility is that the 
effect results from non-uniform system transmission affecting the sodium as 
well as the background radiation. Indeed, as noted below, many of our back- 
ground-subtracted images show generally brighter emission to the north of lo 
along a north-south line passing through the satellite. However, others show 
brighter emission to the south and some show no appreciable asymmetry at all 
along this line (e.g.. Figure 7d). We conclude that the north-south variations 
in the brightness of the sodium cloud at lo's distance from Jupiter in our data 
reflect primarily previously observed effects (e.g., Trafton, 1977, 1980), dis- 
cussed below, of the oscillating plasma torus. 

We made no attempt to observe photometric standards during these observa- 
tions. However, it is possible to obtain a rough estimate of the absolute 
brightness levels of data that fell within the reproducible range of our inten- 
sity calibration curves. We assume that the brightest [S II] X6731 emission in 
the images discussed by Pilcher et al. (1981) is 800 R, a value in accord with 
photometric observations of this line (Trauger et al . , 1980; Brown and Sheman- 
sky, 1982). Inspection of contour plots of the brightest of these data, cor- 
rected for the system spatial transmission function, yields a corrected- 
relative-intensity (CRI)-to-Raylelghs conversion factor of ~1 . 1 . Allowance for 
the fact that the sulfur exposures were 10 minutes while those for sodium were 
only 3 minutes leads to a CRI-to-Raylelghs conversion factor of for the 
sodium data. (The conversion factor for uncorrected data varies across the 
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field of view.) It should be kept in mind, however, that for uncorrected rela- 
tive intensity values above ~400 (CRI values above ~1600), the intensity scale 
may be substantially non-linear. Accordingly, CRI values much above ~1600 
should be regarded as having little quantitative significance. 
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Ill . OBSERVATIONAL RESULTS 


The observations discussed here were obtained in 1980 and 1981. The dates 
and times are given in Table I. The aspect of the data we will address in this 
paper is the presence of directional features, defined as elongated features 
whose symmetry axes are generally inclined to the direction expected for fea- 
tures in lo’s orbital plane. We will further confine our discussion to fea- 
tures directed away from Jupiter as projected onto the plane of the sky. Since 
these features usually extend beyond lo's orbital distance of 5.9 Rj, they must 
comprise sodium whose initial speed at 5.9 Rj is greater than lo's orbital 
speed. Features directed Inward toward Jupiter, also present in our data, may 
result from ejection in other directions and at lower velocities. 

Outward directional features, inclined both to the north and south of lo's 
orbital plane, were detected on almost every occasion on which we observed in 
1980, but less often in 1981. On one occasion, January 6 UT, 1980, a feature 
was seen to change its direction from southward to northward. Three images of 
this feature are illustrated in Figures 3 and 4. (In these observations the 
variation in system transmission across the field-of-view resulted in severe 
attenuation of sodium emission Inward of lo.) In the first observation shown 
in Figures 3a and 4a, the directional feature is inclined slightly toward the 
south. Less than two hours later (Figures 3b and 4b), the feature was inclined 
slightly toward the north. Approximately another two hours later (Figures 3c 
and 4c), the observed feature's axis of symmetry was inclined even farther to 
the north. A background image obtained as part of this sequence is shown in 
Figure 3d. Any instrumental effect associated with the weak north-south inten- 
sity gradient present in the background cannot be responsible for the much 
stronger north-south variations manifested by the directional features. Other 
images obtained in this sequence showed the sodium emission inward of lo's 
position to be concentrated to the north of lo's orbital plane throughout the 
~ 4-hour Interval of observations. The directionality of the final northward 
feature shown in Figures 3c and 4c may therefore be somewhat exagerated, since 
an image showing all of the sodium emission present at that time would likely 
show a broad region of brighter emission to the north. 

The general morphology of the feature observed in this sequence, particu- 
larly as seen in Figures 3a and 4a, is not unusual. In Figure 5 we show a slm- 



ilar feature, although inclined to the north, observed more than a year later. 
We therefore proceed under the assumption that the change in inclination of the 
apparent symmetry axis illustrated in Figures 3 and 4 is real and characteris- 
tic in general of the process forming the directional features. This change in 
inclination, and the more general observation that features may be directed to 
either side of lo's orbital plane, indicate that the features are not associ- 
ated with any single region on lo's solid surface. Rather, a process involving 
interaction with the plasma torus, which oscillates north-south with respect to 
lo , is suggested. We discuss such a model in Section IV. 

A large Increase in the brightness of the sodium cloud around lo and in 
the general region of the directional features was observed in February 1980. 

A group of images from the second night of the four-day interval in that month 
during which we acquired data was shown in Figure 1. In Figure 6 we show 
transmission-corrected contour plots corresponding to three of those images. 
Although of much higher intensity, the directional features observed during 
this period share morphological characteristics with some of the rest of the 
data. The directional feature in Figures 1 and 6 is toward the north and 
exhibits a narrow core embedded in more diffuse emission. This latter charac- 
teristic is also shown weakly, and at much lower brightness, in the image il- 
lustrated in Figures 3c and 4c. The emission on the Jupiter side of lo is 
centered on the projection of the leading portion of lo's orbit onto the plane 
of the sky (cf. Fig. 6). This is generally consistent with the description of 
the region-B cloud developed in previous imaging studies (Murcray and Goody, 
1978; Goldberg et al. 1980; cf. review by Pilcher and Strobel, 1982). There is 
also an overall asymmetry to the sodium cloud along a north/south line passing 
through lo, as discussed above, that is best seen in the contour maps of Figure 
6. In both this case and that illustrated in Figure 4, higher emission inten- 
sities are present to the north. Since all of these data were acquired while 
lo was at northerly magnetic latitudes, this asymmetry is likely an effect of 
electron impact ionization of sodium in the plasma torus which was predominant- 
ly to the south of lo. This asymmetry, which was first observed and discussed 
by Trafton and Macy (1975; Trafton, 1977, 1980), has since been reported by 
other observers (Munch and Bergstralh 1977; Murcray and Goody 1978; Pilcher and 
Schempp 19 79). 
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In Figure 7 we show three images acquired on the fourth night of the 
period of enhanced sodium emission around lo. The feature in this case is much 
narrower and directed away from Jupiter toward the south at an inclination on 
the plane of the sky roughly equal and opposite to that of Figures 1 and 6. In 
Figure 7, however, the feature appears to be displaced: its extension across 

the occulting disk does not cross the disk's center, i.e., lo's position. The 
displacement at lo's radial distance from Jupiter is ~5 arcseconds north. The 
feature in Figure 7 is also distinctly non-linear. The inclination of the 
region B cloud, shown projected toward Jupiter most clearly in the contour plot 
of Figure 7d, does not follow the leading portion of lo's orbit as it does in 
Figure 6. Explanations for this inclination, as well as an analysis of other 
information on region B contained in our data, are beyond the desired scope of 
this paper. Two additional examples of southward directional features are 
shown in Figure 8. 

The evidence that the sodium emission intensity was significantly enhanced 
for the entire four day period that data were obtained in February 1980 cannot 
be attributed to geometrical effects such as variations in lo's orbital phase 
angle. Nor does it appear to be due to photographic or instrumental effects, 
since the photographic densities of the background images obtained during 
February are similar to those of background images from January and from later 
observations at comparable zenocentric distances. Our observation of only one 
period of enhanced emission, and the fact that in other studies of essentially 
the same region of circura-Io space no such period has been observed (Murcray 
and Goody, 1978; Goldberg et al. 1980), indicate that the enhancement is a 
relatively infrequent occurence. The absence of enhanced emission in the data 
of Bergstralh et al. (1975, 1977) may be further indication. It may, however, 
reflect the fact that the small lo-centered region observed by Bergstralh et 
al. (which is completely obscured by our occulting disc) is dominated by sodium 
ejected from lo at velocities lower than those necessary for directional fea- 
ture formation. These lower velocity components may be less affected by the 
processes responsible for enhancing the high velocity sodium seen in the direc- 
tional features. A final point to note here is that images of the plasma torus 
obtained in the light of [SII] X6716 on the same nights in February that the 
enhanced directional features were observed showed no unusual characteristics, 
indicating that a change in plasma properties was probably not responsible for 
the increased intensity. If this reasoning is correct, then the rate of sodium 
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escape from lo during this period at the moderate to high velocities required 
to populate the observed cloud (see Section IV) must have been anomalously 
high. 

As noted above, observations of both northward and southward directional 
features, and the one case of a feature that changed its orientation from south 
to north, suggest that a process involving the plasma torus, which oscillates 
about the satellite plane, is involved in the formation of the directional fea- 
tures. To test this suggestion we divided all of our observations into three 
categories: those showing a northward directional feature, those showing a 

southward feature, and those not showing a directional feature. We show the 
results of this categorization in Figure 9 plotted with respect to the magnetic 
longitude of lo at the time of observation. It can be seen that all of the 
northward directional features were observed when lo was between 170 and 360 
degrees magnetic longitude. If the observation of one very weak northward fea- 
ture which faded during the observation is excluded (the data of April 13 UT, 
1981), then all the remaining northward features were observed at lo magnetic 
longitudes between 170 and 300 degrees. Southward features were almost entire- 
ly restricted to longitudes outside this range. On two occasions features dis- 
appeared during an observing period, and on one occasion, discussed at length 
above, a feature changed its direction. The data in Figure 9 indicate that the 
north-south orientation of the directional features is correlated with the mag- 
netic longitude of lo, and that southward features seem to fade or change to 
northward features for lo magnetic longitudes of about 150 to 160 degrees. 
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IV. ANALYSIS 


We will confine our analysis here to an interpretation of the characteris- 
tics of the directional features, defined at the beginning of Section III as 
those features both inclined to the direction expected for features in lo's 
orbital plane and projected onto the plane of the sky in the direction away 
from Jupiter. The bright reglon-B sodium cloud and the occasional peculiar 
sodium structural features visible inside of lo’s orbit in some of the observa- 
tions will not be analyzed here. 

The Implication of the variation with magnetic longitude shown in Figure 9 
is that the features would appear to a terrestrial observer to oscillate north 
and south of the satellite orbital plane with a period of about 13 hours, the 
orbital period of lo relative to Jupiter's magnetic field. This inferred 
oscillation indicates that the features are not produced by some active local 
sodium source on the satellite’s surface. It seems more likely that the oscil-j 
lation results from the modulation of the sodium cloud density near the satel- 
lite plane introduced by the time-dependent lifetime of sodium in the oscillat- 
ing plasma torus. This is the main idea to be pursued in this analysis. 

The lifetime of a sodium atom in the lo plasma torus is determined primar- 
ily by the rate of impact ionization by plasma torus electrons. The spatial 
variation of this lifetime is shown in Figure 10 (Smyth, 1983). Along the ver- 
tical oscillation of lo, also shown in Figure 10, the lifetime of sodium is 
small and varies between about one and four hours. Because of this oscillation 
it is possible for sodium atoms ejected from lo with a significant northward 
velocity component to traverse a region of quite different sodium lifetime than 
atoms ejected simultaneously to the south. An atom ejected to the north when 
lo is at the northernmost extreme of its oscillation, for example, will enter a 
region of comparatively long lifetime, while one ejected simultaneously to the 
south is likely to be ionized relatively quickly. This difference will lead to 
a significant north-south time-dependent asymmetry in the sodium cloud density 
outside of lo’s orbit (the region of interest for the directional features) if 
the atoms have speeds within a certain range so as to move in proper phase with 
the torus’ oscillations. If the speeds are too small, the atoms will require 
half an oscillation period or more to move a significant distance in the plasma 
torus, and the lifetime histories for the north and south sodium will have 
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similar time-averaged values. If the speeds are too large, the north and south 
sodium will again have similar lifetime histories since neither will spend 
enough time in the torus to be significantly affected. If the directional fea- 
tures are to be explained in terms of the oscillating sink, the relevant sodium 
source must involve only atoms having an intermediate range of initial speeds. 

To define better this intermediate range, one can use the criterion that 
the initial velocities of sodium atoms relative to the satellite must be large 
enough to at least populate the circular viewing window in a time interval ap- 
proximately one-quarter as long as the 13 hour oscillation period of the plasma 
torus with respect to lo. Several relevant orbit trajectories are shown in 
Figure 11 projected upon the sky plane for the observational geometry of Plate 
1076a (see Figures 3c and 4c). In both Figures Ha and 11b, the three dif- 
ferent trajectories correspond to initial velocity vectors that are inclined 30 
degrees above the satellite plane and have initial speeds relative to lo of 6, 
10 and 20 km sec~^. The components of the initial velocity vectors in the 
satellite plane are directed along the orbital motion of lo for Figure 11a and 
along the anti-Jupiter direction for Figure 11b. These trajectories suggest 
that initial speeds of approximately 20 km sec”^ satisfy the speed criterion, 
with slightly larger speeds acceptable for the initial velocity direction of 
Figure 11a and slightly smaller speeds acceptable for that of Figure 11b. Ini- 
tial speeds of 10 km sec ^ or less are too slow and more than 30 km sec ^ (not 
shown in Figure 11) are too fast. 

With the speed range defined, we used the lo sodium cloud model of Smyth 
(1983), recently improved to include the oscillating plasma torus sink for 
sodium (Smyth and Combi, 1983; also see Smyth, 1983), to perform several ex- 
ploratory calculations. The observational sequence shown in Figures 3 and 4 
was selected for a model-data comparison, since it provides the potentially 
diagnostic observation of a change in feature direction from south to north 
(these observations are also illustrated as the January 6 UT, 1980 data in 
Figure 9). In our first exploratory model calculations, we assumed the sodium 
to be ejected at 20 km sec ^ from the trailing hemisphere exobase of the satel- 
lite (excluding a stagnation region centered on the trailing hemisphere apex). 
Our choice of the trailing hemisphere to be the source region was motivated by 
the thought that the directional features may result from an interaction 
between lo and the corotating magnetospheric plasma which most directly im- 
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pinges on the satellite's trailing hemisphere. At the assumed velocity of 
20 km sec solar radiation pressure effects can be ignored. The initial 
velocity vectors on the exobase were not radially directed, but were arranged 
in a cylindrically symmetric fashion about the direction of lo's orbital motion 
and inclined by an angle a above the forward-directed tangent vector to the 
exobase surface. The angle a was assumed to have small values, mostly between 
0 and 20 degrees, as suggested by line profile studies (cf. Section I) which 
showed sodium to be ejected in the forward direction and by the momentum direc- 
tion of the raagnetospheric plasma moving past lo. 

The exploratory model results for the three observations of Figures 3 and 
4 did not produce distinct directional features. Furthermore, the north-south 
asymmetry in the calculated sodium Intensity on the sky plane was toward the 
south for both Plates 1069c and 1073a and only weakly toward the north for 
Plate 1076a. This is in disagreement not only with this particular sequence of 
observations, but also with the general trend of feature directionality with lo 
magnetic longitude illustrated in Figure 9. It can be seen from this figure 
and, for example, Figures 1 and 6, that Plate 1076a was acquired at an lo mag- 
netic longitude well into the range in which strongly northward directed fea- 
tures have been observed. The quite weak northward enhancement calculated for 
this plate thus seems to be generally inconsistent with the data. The situa- 
tion was not significantly improved when we considered a severe system III mag- 
netic longitude dependence to the electron impact ionization lifetime of sodium 
in the plasma torus, nor when we included an initial speed dispersion for the 
sodium atoms. This suggests that the north-south orientations of the direc- 
tional features are controlled more by the general oscillatory nature of the 
plasma torus than by the fine details of the spatial lifetime structure of 
sodium. It also suggests that the orientations of the initial velocity vectors 
chosen for the exploratory calculations are not appropriate to explain the di- 
rectional features. 

In order to search more effectively for the proper orientations of the 
initial velocity vectors for the sodium atoms, we developed a correlation tech- 
nique and applied it to the observations of Figures 3 and 4. Our aim was to 
correlate the initial velocity directions of sodium atoms with their resulting 
contributions to the D-line intensity in the observing window. To cover all 
possible directions, it was sufficient to consider a set of initial velocity 



vectors that were directed radially and Isotropically from lo's exobase. This 
approach is valid for the relatively large initial speeds of 20 km sec~^ under 
consideration, since the path of an atom leaving lo at such a speed is almost 
totally dependent upon its initial velocity direction and almost totally inde- 
pendent of its initial location on the exobase. For purposes of this study, 
the exobase of lo (assumed to be 2600 km in radius) was divided into 613 north- 
south pairs of points, each pair comprising two equal areas arranged symmetri- 
cally with respect to the equator. The central latitudes and longitude of each 
pair were taken to be the coordinates of radial ejection of sodium atoms from 
the exobase. The lo sodium cloud model was then used to calculate the space- 
time D 2 -line intensity contribution along each north-south pair of trajectories 
for an initial speed of 20 km sec ^ and the magnetic and viewing geometries 
corresponding to the observations of Figures 3 and 4. 

For the purpose of this study, we defined upper and lower windows in the 
main observing window as shown in Figure 12. For each trajectory pair, we cal- 
culated the cumulative D 2 -Hne intensity contributions to the upper window ly 
and the lower window ly. We found that a 20-hour flight time was the maximum 
we need consider. We also computed the Intensity ratio given by 


R 



( 1 ) 


If R is greater than one, the trajectory pair contributes primarily to a 
northern enhancement in the D-line intensity. If R is less than one, the en- 
hancement is to the south. For each trajectory pair, the value of R can there- 
fore be associated with two angles, the ejection longitude and the absolute 
ejection latitude. Since the atoms are assumed to be ejected radially, these 
angles describe the initial velocity directions of the pair. Our results for 
all of the 613 north-south pairs that make a contribution to the intensities in 
the upper and lower viewing windows are summarized in Figure 13. We have fol- 
lowed here the normal longitude convention in which the sub-Jupiter point is at 
0° and the apex of motion is at 90°. The three parts of Figure 13 correspond 
to the geometries of the three observations of Figures 3 and 4. Values of R 
greater than unity have been displayed as three sizes of the letter N, while 
corresponding values of R ^ have been displayed as the same three sizes of the 
letter S. The large N (or S) symbols correspond to values of R (or R larger 

than 2, the intermediate N (or S) symbols to values between 2.0 and 1.5, and 
the small N (or S) symbols to values between 1.5 and 1.0. 
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The results of the correlation study in Figure 13 are striking. These 
results show that sodium atoms with initial velocity vectors directed radially 
outward and somewhat normal to lo's orbital motion (i.e., 180° ejection longi- 
tude) have the correct orientation to explain generally the south-to-north 
changes that occur in the directional features shown in Figures 3 and 4. Re- 
sults in Figure 13 also confirm our earlier exploratory model conclusions that 
sodium atoms ejected with initial velocity vectors directed nearly parallel to 
lo’s orbital motion produce intensity patterns that are too bright in the south 
(see the large S's at low latitudes and longitudes less than 90° in Figure 13). 

The spatial morphology of the sodium D 2 -llne intensity seen on the sky 
plane through the viewing window, produced by 20 km sec ^ sodium atoms directed 
at nearly right angles to lo's orbital motion, is therefore the appropriate 
situation to be explored. To accomplish this we used the lo sodium cloud model 
with the ejection of sodium atoms from lo restricted to the band illustrated in 
Figure 14. The band is centered on the great circle that is normal to lo's 
orbital motion direction and extends about 20 degrees toward the leading and 
trailing hemisphere apexes. In the model calculation 456 trajectories, corres- 
ponding to atoms emitted from equal areas distributed symmetrically around the 
band, were used to simulate the escape of sodium from lo. For equal sodium 
flux from all areas and a value for the ejection angle a of 90 degrees (i.e., 
radial ejection), the calculated D 2 intensity distributions corresponding to 
the observational geometries of Figures 3 and 4 are shown in Figure 15. Note 
that in Figure 15 the general intensity pattern calculated from the model 
changes from slightly south to north in phase with the observations. The model 
calculation thus exhibits the correct dependence on the magnetic longitude of 
lo, but lacks the distinct spatial character of the observed features. 

The extreme northern enhancement of the calculated sodium intensity for 
Plate 1076a in Figure 15 is to be expected and may be understood as follows. 

As discussed above, sodium atoms ejected from lo to the north at the intermedi- 
ate velocities in question have a decided lifetime advantage if they are 
ejected when lo is near the northern extreme of its oscillation with respect to 
the plasma (cf. Fig. 10). A northern enhancement in the sodium cloud density 
then results. This is the case (see Fig. 13c) for Plate 1076a (Figures 3c and 
4c), which was acquired when lo was at 9.2° north magnetic latitude and had 
been at high northern magnetic latitudes for the preceding 2-3 hours. A cor- 
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responding southern enhancement results in the calculations when lo is assumed 
to be near the southern extreme of its oscillation. The overall effect in the 
model is a broadening of the emission distribution to the north or south with a 
corresponding loss of specific directionality. 

To obtain more pronounced and definitely directional features in the 
model, as are seen in the observations, this high latitude lifetime advantage 
must be compensated for by an enhancement of the escape flux of sodium from low 
latitudes. We therefore examined the consequences of a latitude-weighted dis- 
tribution in which the flux of sodium (j) from the band ejection region is des- 
cribed by the expression 

<() = (}) (1 - f |sin ijj ), (2) 

0 

where (j) is the flux on the band at the equator, $ is the angle measured from 
0 

the equator along a line parallel to the great circle that is the center of the 
band ejection region, and f is a reduction factor chosen from a model-data com- 
parison to have a value of 0.7. 

The results from this latitude-weighted band ejection model, calculated 
for the observing geometry of Plate 1069c (Figures 3a and 4a), are shown in 
Figure 16 for several values of the ejection angle a and an initial speed of 
20 km sec”^. The results for ejection angles of 75°, 90°, and 105° show a di- 
rectional character to the emission outward of lo that is more consistent with 
the observations than the unweighted results. The intensity distribution cal- 
culated for 105° most clearly shows southward directionality, in accord with 
the corresponding observation. 

For ejection angles of 60 degrees and smaller, a more prominent southward 
directional feature forms that exhibits neither the correct spatial morphology 
nor the observed lo magnetic longitude phase relation shown in Figure 9. This 
is consistent with our earlier exploratory runs and with the correlation 
results of Figure 13. For ejection angles of 120° and larger, it becomes in- 
creasingly difficult to produce significant directionality. This occurs for 
two reasons. First, the lower resultant orbital velocities of atoms ejected on 
these trajectories result in their not traveling as far outside of lo’s orbit 
as atoms ejected at smaller angles, and therefore not being seen as far from 
lo's position on the plane of the sky. Second, at these smaller distances from 
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Jupiter, the atoms are more readily ionized by the plasma and hence contribute 
less D-line emission intensity. It therefore appears that the directions of 
the initial atomic velocity vectors most consistent with the morphologies of 
the directional features are outward from lo, at approximately right angles to 
the orbital motion of the satellite, with atomic flux generally decreasing from 
equator to pole. 

The results of model calculations made with this weighted band ejection 
model for all three of the observations of Figures 3 and 4 are shown in Figure 

17. Here we have adopted an ejection angle a of 90° (i.e., radial ejection), 
although the use of 75° or 105° gives very similar results. The calculated 
intensity distribution outward of lo is similar in general shape and dependence 
on lo magnetic longitude to the observations. 

Since the directional features in these model calculations result from 
initial velocity vectors directed outward from lo's orbit within a narrow lon- 
gitudinal angular range, the question arises as to whether the features could 
just as well be produced by radial ejection from the trailing hemisphere of 
lo. This may be possible since radial ejection from the portion of the trail- 
ing hemisphere that does not overlap the band region is less able to contribute 
to intensity within the observing window. There are two reasons for this. 
First, the atoms on the trajectories that originate on the trailing hemisphere 
outside the band region lie closer to lo's orbit and hence are more susceptible 
to ionization by the plasma torus. Second, a significant fraction of these 
atoms rapidly move inside of lo's orbit and are either lost from the cloud by 
collision with Jupiter or contribute only to the near-Io region-B cloud. To 
pursue this idea, we performed calculations in which the trailing hemisphere 
was divided into 649 equal area segments, from each of which we assumed sodium 
to be ejected radially at 20 km sec”^. The escape flux ij) of sodium was en- 
hanced near the equator through the use of a latitudinal distribution given by 

expression (2) where ({) is the sodium flux at the equator, $ is now latitude, 

0 

and f is a reduction factor chosen again to have a value of 0.7. 

Results for all three observations of Figures 3 and 4 are shown in Figure 

18. The general spatial morphology of these calculated D 2 intensity contours 
compares favorably with the observations. For Plate 1069c, however, the calcu- 
lated feature is not as vertically distinct as in the observation (Figure 4) or 
as in the band ejection model results (Figure 17). This reduced directionality 


21 



is apparent also in the band ejection results of Figure 16 for an ejection 
angle a of 120 degrees, and reflects the importance of velocity vectors on the 
forward half of the band ejection region. Any real ejection region centered on 
the trailing hemisphere of lo would, of course, likely yield some velocity vec- 
tors with forward components. This would naturally extend the effective exo- 
base region of ejection beyond the strict geometric boundary adopted in the 
model. As anticipated, apart from these differences, the spatial intensity 
morphologies of the features shown in Figure 18 are similar to those calculated 
for the band ejection and shown in Figure 17. 

From examination of the correlation results of Figure 13, one might also 
ask if the directional features of Figures 3 and 4 could be suitably modeled by 
radial ejection from the leading hemisphere. This may be possible if ejection 
from that portion of the leading hemisphere that overlaps half of the band re- 
gion is sufficiently strong to dilute the intensity contributions of the ejec- 
tion originating closer to the hemispherical apex, which has been shown earlier 
to provide a southern enhancement. To examine this situation, we simulated 
leading hemisphere ejection from the satellite exobase by a procedure corres- 
ponding exactly to that adopted for the trailing hemisphere. The model results 
for all three observations of Figures 3 and 4 are shown in Figure 19. The 
general change in the direction of the feature from south to north is in the 

correct sense but can be seen to lag behind the observation (i.e., the direc- 

tion of the calculated feature advances toward the north too slowly). The 
north-south character of the feature is also broader than in the observations 
or the band ejection results (Figure 17). These characteristics result as an- 
ticipated, from the southern intensity contributions made by sodium ejected 
near the apex of the leading hemisphere and around the leading-inner quadrant 
at low latitudes (see Figure 13). An appropriate variation of atomic flux with 
longitude that reduced the importance of these atoms could therefore be used to 
improve the agreement of this variation of the model with the observations. 
However, the presence of some of these near forward-directed atoms with speeds 

up to 20 km sec”^ has been shown to be required for the proper modeling of the 

asymmetric skirt of the sodium D-line profile (see the discussion in Section 
I). In a model that accounts for both the directional features and the D-line 
profiles, therefore, the directional features may prove to be understandable as 
the spatial counterpart of the spectral signature given by the asymmetric 
skirt . 
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By comparing the estimated brightness levels from Figure 4 with the model 

results of Figures 17, 18 and 19, we may estimate the satellite equatorial 

atomic flux (J> and the total rate of sodium ejection from lo at these veloci- 
0 

ties. As discussed in Section II, the brightnesses corresponding to the con- 
tours plotted in Figure 4 may be estimated (within a factor of 2-3) by means of 
the rough conversion factor of ~4 from CRI (corrected relative intensity) units 
to Rayleighs. For the band ejection models of Figure 17, the equatorial flux (j> 
required to match the observations is about 2.5 x 10^ atoms cm“^ sec”^ for 
Plate 1069c and about 2.1 x 10^ atoms cm~^ sec~^ for Plates 1073a and 1076a. 
Integrating the non-uniform flux distribution over the band region, we find 
this corresponds to a total satellite source rate between about 1.7 and 2.0 x 
10 ° atoms sec~^. Only about half of these atoms contribute to the directional 
features in the viewing window (see Figure 13). For the trailing hemisphere 
ejection model results of Figure 18, the required equatorial flux is 5 x 10^, 
3.2 X 10^ and 2.7 x 10^ atoms cm~^ sec~^, respectively, for Plates 1069c, 

1073a, and 1076a. This corresponds to a total satellite source rate between 
about 3.7 and 6.8 x 10^^ atoms sec“^ for the postulated non-uniform flux dis- 
tribution over the trailing hemisphere exobase. At most half of this material 
contributes to the brightness of the directional feature seen in the viewing 
window. For the leading hemisphere ejection model results of Figure 19, the 

Q —9 —1 ’ 

required equatorial flux is about 1.4 x 10 atoms cm sec for all three 
observations. This corresponds to a total satellite source rate of about 1.8 x 
10 ° atoms sec ^ of which slightly more than half may contribute to the direc- 
tional feature. From these comparisons, we conclude that the source rate 
required to produce the directional features is of order 10“^^ atoms sec 
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V. DISCUSSION 


The general conclusions that we have derived in the analysis may be stated 
as follows. The directional features result from the combined effects of a 
source of high speed (~20 km sec~^) sodium and the oscillating lo-plasma-torus 
sink for neutral sodium. The lo magnetic longitude at which the computed fea- 
tures in our model change direction with respect to lo's orbital plane is cri- 
tically dependent on the directions of the initial velocity vectors of the high 
speed atoms. To account for the observations, these must be directed essen- 
tially at right angles to lo's orbital motion. The directionality of the fea- 
tures requires that the flux of sodium escaping from lo generally decrease from 
equator to pole. The required initial velocity directions are consistent with 

(1) radial atomic ejection from a band region of narrow longitudinal extent 
centered on the great circle separating the leading and trailing hemispheres, 

(2) approximately radial ejection from the trailing hemisphere (possibly 
extending forward into the leading hemisphere by 20° and possibly with little 
or no ejection from a region centered on the trailing hemisphere apex), or (3) 
radial ejection from the leading hemisphere with a substantial decrease in flux 
toward the apex of motion. The actual sodium source may be distributed over 
some combination of these regions, as illustrated in Figure 20a. 

We have not thus far appealed to any particular physical mechanism to ex- 
plain the deduced characteristics of the sodium ejection from lo. We note, 
however, that the required right-angle direction of ejection of sodium from the 
satellite relative to its orbital motion is exactly the same direction shown to 
occur for the elastic scattering of neutral atoms from lo by magnetospheric ion 
impacts (Sieveka, 1983; Sieveka and Johnson, 1984). The near right angle ori- 
entation of the neutral atoms follows naturally from the scattering properties 
of a screened coulomb potential which realistically describes the ion-neutral 
collisional interaction. This is in marked contrast to the highly forward- 
peaked scattering which results from hard-sphere interactions. The required 
single scattering of atoms could occur in an extended neutral corona of the 
satellite, producing a substantial atomic flux in the 10 to 30 km sec~^ velo- 
city range. The corona may be populated by ion sputtering of neutrals from the 
trailing hemisphere of lo's surface (in the low density atmospheric limit) or 
exobase (in the high-density atmospheric limit). In addition, somewhat radial 
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ejection of neutral atoms may be produced directly by sputtering from the 
trailing hemisphere surface or exobase. 

Our conclusions, based on analysis of the directional features together 
with the calculations of Sieveka (1983) and Sieveka and Johnson (1984), provide 
direct evidence for a magnetospheric-wind-driven gas escape mechanism. Either 
the somewhat radially directed sodium atoms expected from sputtering, or the 
near right angle source of sodium provided by scattering in the corona (or, 
realistically, a combination of the two) may supply the necessary atoms to ex- 
plain the directional features. The relative importance of the sputtered and 
corona sources will depend upon the local properties of the satellite atmos- 
phere (and/or surface) and its interaction with the plasma flow pattern in the 
vicinity of lo. A simple flow pattern that is suggested by the ejection pat- 
tern of Figure 20a is illustrated in Figure 20b. More elaborate plasma flow 
patterns based upon a two-dimensional electric dipole have been adopted by Ip 
(1982) and used to argue that such a dynamic coupling of a neutral atmosphere 
and the corotating plasma can provide a significant source for the neutral 
clouds in the Jovian system. 

Further quantitative study of the directional feature data coupled with 
calculations of the flux and angular distribution of escaping material from the 
satellite surface or exobase and from its corona may allow definite conclusions 
to be drawn as to the presence or absence of a substantial atmosphere on lo. 

In addition to this information, the indication that the sodium flux is en- 
hanced at the equator may reflect a latitudinal structure in the atmosphere, 
perhaps produced by volcanic activity near the equator or by colder tempera- 
tures nearer the satellite poles. The very strong directionality of some fea- 
tures not analyzed here (see, for example. Figures 1, 6, and 7) may reflect 
temporal variation in this structure. Observations in which the directional 
features are weak or absent may also be diagnostic of changes in the magneto- 
spheric properties or satellite atmosphere. Further studies of the directional 
feature data together with the region-B cloud data (not considered here) are 
necessary to determine the total source rate of sodium. This is important for 
understanding the source rates and escape mechanisms operative at lo for other 
gases such as those containing oxygen and sulfur, the ions of which constitute 
the bulk of the plasma in Jupiter's inner magnetosphere. 
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Although the direct evidence deduced from the directional feature data 
supporting a magnetospheric-wind-driven gas escape mechanism for lo is new, the 
general subject has been under discussion for several years. Charged-particle 
sputtering of sodium from lo's surface was first suggested by Matson et al. 
(1974). Laboratory investigations of ion sputtering from various icy surfaces 
have subsequently been actively pursued (Brown et al. 1980; Lanzerotti et al. 
1982; Brown, Lanzerotti and Johnson, 1982; Johnson et al. 1983) to evaluate the 
mechanism more quantitatively. Sputtering from the atmosphere of lo by the 
corotating thermal plasma was also initially discussed by Haff, Watson and Yung 
(1981) and more recently by Summers, Yung, and Haff (1983). Single-scatter 
elastic collisions of neutrals with the corotating thermal plasma (i.e., knock- 
on processes) and charge exchange processes have been discussed by several 
authors (Brown and Schneider, 1981; Ip, 1982; Brown, Pilcher and Strobel, 

1983). The recent efforts of Sieveka (1983) and Sieveka and Johnson (1984) 
have been more specific and have begun to address these various escape 
processes collectively and more quantitatively. Additional efforts in this 
direction should prove to be very valuable in our attempts to understand better 
the characteristics of the directional features of the sodium cloud. 
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Table I 


Date(UT) 


Plate Numbers 


Start Time 
(UT) 


Jan. 5 

1065,1066 

1207 

Jan. 6 

1069-1076 

1014 

Jan. 7 

1079 

1144 

Feb . 5 

1083 

842 

Feb . 6 

1088-1091 

823 

Feb. 7 

1095-1098 

946 

Feb. 8 

1103-1105 

836 

Apr. 21 

1191-1194 

922 

Apr . 22 

1203 

921 

Mar. 10 

1446-1454 

722 

Mar . 11 

1455,1456 

821 

Mar . 1 1 

1459-1462 

1101 

Mar . 12 

1465,1466 

833 

Mar. 12 

1471,1472 

1351 

Mar. 13 

1473,1474 

700 

Apr. 9 

1838-1841 

1016 

Apr. 10 

1845,1846 

909 

Apr. 13 

1862-1869 

551 


End Time 
(UT) 


1303 

1424 

1153 

852 

1000 

1119 

947 

1050 

933 

1439 

850 

1323 

911 

1420 

730 

1203 

1021 

1004 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Photographs of plates of the sodium emission around lo obtained on 
1980 Feb. 6 between 823 and 909 UT. Each image is 110 arcseconds 
(5.0 Rj) in diameter. lo's image is centered on the lA-arcsecond 
occulting disk that appears as a bright spot in the center of each 
negative image. North is up and east is to the left. Jupiter is 
to the left, lo being essentially at western elongation when these 
images were acquired (see Fig. 6). The sequence ^ - f_ corresponds 
to a sequence of increasing filter tilt angle from 1° to 6° in 1° 
increments. The differences between images are explained in the 
text. The vertical intensity discontinuity near the left edge of 
each image, seen best in images c_ - J_, is an artifact caused by 
scattering within the telescope. The directional feature, evident 
in images ^ extends from the occulting disk toward the upper 

right to the margin of the field. Image is a measurement of the 
background which was used in the reduction of images ^ (see 
text ) . 

(a) The calculated instrumental spatial transmission function for 
sodium D 2 ~line radiation corresponding to the image of Fig. lb, 
which was acquired at a filter tilt angle of 2°; (b) Same as but 
for the image of Fig. Ic, which was acquired at a filter tilt angle 
of 3°; (c) Same as a_, but for the image of Fig. Id, which was 
acquired at a filter tile angle of 4°. 

Photographs of plates obtained on 1980 Jan. 6 UT. North is up and 
east is to the left. Jupiter is to the right. Additional details 
are given in the caption to Fig. 1. Details regarding images - 
which show the directional feature, are given in the legends of 
Fig. 4. In these images, all of which were obtained at a filter 
tilt angle of 2-1/2°, the sodium emission inward of lo (to the 
right of the occulting disk) was severely attenuated owing to the 
spatially non-uniform system transmission. Image ^ is a background 
image obtained 12 minutes after image ^ at a filter tilt angle of 
8 °. 

Contour plots of the images shown in Figures 3a - c. North is up 
and east is to the left. Jupiter, lo's orbit, and lo's position at 
the time of the observation are illustrated. The background has 
been subtracted and the resulting images corrected for the system 
spatial transmission function as described in the text. Inward of 
lo's position the spatial transmission function for these 
observations is so low that little sodium radiation is 
transmitted. The contour values in corrected-relative-intensity 
(CRI) units are as follows: (a) 65, 112, 150, 250, 450, 750, 1100; 

(b) 65, 90, 200, 500, 1200; (c) 65, 100, 190, 400, 800, 1400. For 
a rough conversion from CRI units to Rayleighs, multiply by a 
factor of 4. 

Photographs of images obtained on 1981 April 10 between 909 and 916 
UT. Jupiter is to the right (west). See the caption to Fig. 1 for 
additional details. The unocculted image of Ganymede, to the lower 
right of the occulting disk, is surrounded by rings of scattered 
light, one of which transects the directional feature. 
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Figure 


Figure 


Figure 


Figure 


Figure 


6. Contour plots of the images shown in Figures lb - d. The format 
and data processing are as described in the caption to Fig. 4. The 
contour values in CRI units are as follows: (a) 65, 150, 250, 375, 

625, 1000, 1750, 3500, 5000, 6500; (b) 65, 150, 250, 375, 625, 

1000, 1750, 3500, 5000; (c) 65, 150, 250, 375, 625, 1000, 1750, 

3500, 5000. As discussed in the text, CRI values in excess of 
~1600 should be regarded as having only qualitative significance. 

The CRI values of these three images are in good agreement despite 
substantial differences between the uncorrected intensity levels. 

7. (a-c) Photograph of a plate obtained on 1980 Feb. 8 between 836 and 
850 UT. Jupiter is to the left (east). See the caption to Fig. 1 
for additional details; (d) A contour plot of the image shown in 

The format and data processing are as described in the caption 
to Fig. 4. The contour values in CRI units are as follows: 65, 

100, 150, 250, 375, 625, 1250, 2500, 5000. CRI values above -1600 
should be regarded as having only qualitative significance. 

Multiplication by a factor of 4 roughly converts lower CRI values 
to Rayleighs. 

8. Photographs of images obtained on 1981 Mar. 11 UT (a^) and 1981 Apr. ^ 

9 UT (_^) showing southward directional features. Jupiter is to the 

right in and to the left in See the caption to Fig. 1 for 
additional details. 

9. Classification of directional features as a function of lo's mag- 

netic longitude at the time of the observations. Observed north- 
ward features are plotted at the top, southward features at the 
bottom, and null observations (no feature) in the center. Obser- 
vations at overlapping longitudes have been displaced vertically 
for clarity. Observations from 1980 are plotted above those from 
1981. Solid lines connect data points corresponding to a con- 
tinuous series of observations on a single night. A dashed line : 

connects the points corresponding to two observations obtained 

—3 hours apart on March 11 UT, 1981. The dotted line partially 
used for the data of April 13 UT, 1981, is meant only to call 
attention to the fact that during these observations lo passed 
through 0° magnetic longitude, and the plotted points therefore 
appear separated in the figure, as do the points for March 11 UT, 

1981. 

10. Sodium Electron Impact Ionization Lifetime in the lo Plasma 
Torus. The lifetime calculation is based upon the cross-sections 
measured by McFarland (1965), McFarland and Kinney (1965) and 
Zapesochnyi and Aleksakhin (1968), Voyager 1 in situ measurements 
of the plasma density (Bridge, Sullivan and Bagenal, 1980; Bagenal 
and Sullivan, 1981), and Voyager 1 UVS (Shemansky, 1980) and in 
situ measurements (Scudder, Sittler and Bridge, 1981) for the elec- 
tron temperature. The two-dimensional common temperature model of 
Bagenal and Sullivan (1981) was selected to describe the electron 
density. 
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Figure 11. Orbit trajectories for Different Initial Velocities. The initial 
speed and total flight time for each trajectory are given, with 
one-hour tick marks indicated. Also shown to scale on the sky 
plane are Jupiter, lo's orbit, lo at 121° orbital phase, and the 
outer boundary of the observing window which is centered on lo. 

See text for orientation of the initial velocity vectors. 

Figure 12. Correlation Study for the Directional Feature. The upper and lower 
segments of the viewing window adopted in the correlation study are 
shown to scale with Jupiter and lo's orbit, for the satellite at 
eastern elongation. 

Figure 13. Results of a Study of the Correlation Between the Angle at which 

Sodium Atoms are Ejected from lo and the North-South Orientation of 
the Resulting Region of Enhanced D-line Emission. Three sizes of N 
or S, corresponding to more northerly or southerly intensity en- 
hancements in the upper (northern) and lower (southern) viewing 
windows defined in Fig. 12 are shown on a latitude-longitude plane. 
The location of a point on this plane corresponds to the location 
on lo's exobase from which the atoms producing the indicated en- 
hancement would be ejected. The ejection of atoms from the exobase 
is assumed to be radial at 20 km sec~^. See the text for more 
quantitative details. 

Figure 14. Band Ejection Model. The directions of the initial velocity vec- 
tors for sodium atoms ejected from the band region, which is 
oriented symmetrically about the direction of lo's orbital motion, 
is illustrated. The angle a can be adjusted to simulate both for- 
ward and backward ejection. 

Figure 15. Results of Unweighted Band Ejection Model Calculations. The calcu- 
lated D 2 ~line intensity contours corresponding to the three obser- 
vations of Figures 3 and 4 are shown. The spacings of the contours 
are the same as those of Figure 4, and the flux of atoms from lo 
has been adjusted to yield approximately the absolute intensities 
estimated from the observations. 

Figure 16. Results of Latitude-Weighted Band Ejection Model Calculations. The 
D 2 “line intensity contours for Plate 1069a (Figure 4a) are shown 
for five values of the ejection angle a. The contour levels were 
determined as described in the caption to Fig. 15. 

Figure 17. Results of Latitude-Weighted Band Ejection Model Calculations. See 
the caption to Fig. 15. 

Figure 18. Results of Latitude-Weighted Trailing Hemisphere Model Calcula- 
tions. See the caption to Fig. 15. 

Figure 19. Results of Latitude-Weighted Leading Hemisphere Model Calculations. 
See the caption to Fig. 15. 
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Figure 20. Magnetospheric-Wind-Driven Gas Escape from lo. (a) The arrows de- 
pict the pattern of atomic ejection vectors from lo inferred in 
this paper to account for the directional features. The most im- 
portant ejection directions, inferred to be the primary source of 
sodium for the directional features, are shown in heavier shading 
for the case of radial atom ejection from lo. The essential 
character of this primary source is its near orthogonal ejection 
direction with respect to the satellite motion, and need not neces- 
sarily be identified directly with a specific region of the exobase 
surface as depicted. (b) A simple plasma flow pattern past lo that 
might produce the ejection pattern in ^ by ion-neutral Interactions 
is shown. 
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